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SUMMARY 

The expression for eddy diffusivity from a previous analysis was 
modified in order to account for the effect of kinematic viscosity in 
reducing the turbulence in the region close to a wall. By using the 
modified expression, good agreement was obtained between predicted and 
experimental results for heat and mass transfer at Prandtl and Schmidt 
numbers between 0.5 and 3000. The effects of length-to-diameter ratio 
and of variable viscosity were also investigated for a wide range of 
Prandtl numbers . 


INTRODUCTION 

Most of the existing analyses for turbulent heat and mass transfer 
are adequate only for Prandtl and Schmidt numbers on the order of 1 or 
less. For instance, the analysis given in reference 1, although adequate 
for gases, gives heat- and mass -transfer coefficients for liquids at high 
Prandtl or Schmidt numbers that are higher than those obtained experimen- 
tally. The difference between the experimental values and the values ob- 
tained by the method in reference 1 increases as the Prandtl or Schmidt 
number increases. Coefficients obtained from the von Karman analysis 
(ref. 2) at high Prandtl or Schmidt numbers are lower than the experi- 
mental values. Rannie's analysis (ref. 3) gives coefficients which are 
in somewhat better agreement with the data than either of these analyses, 
but the coefficients are again inaccurate at very high Prandtl or Schmidt 
numbers. The analysis of reference 4 agrees with data at Prandtl or 
Schmidt numbers of 1 and at very high Prandtl or Schmidt numbers, but the 
coefficients are somewhat low at intermediate values of these numbers . 

In reference 5, which was published since the present investigation was 
initiated, good agreement was obtained with mass -transfer data for 
Prandtl and Schmidt numbers between 0.5 and 3000 by introducing an 
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appropriate amount of turbulence into the laminar sublayer. In all 
these analyses, except those in references 1 and 3, the properties were 
constant and the flow fully developed. The relations among heat trans- 
fer, mass transfer, and fluid friction are discussed in reference 6. 

The inadequacy of most of the previous analyses at high Prandtl and 
Schmidt numbers is principally caused by the expressions used for the 
eddy diffusivity in the region very close to the wall. This region is 
important because of the extremely large temperature or concentration 
gradients in that region at high Prandtl or Schmidt numbers (ref. 7). 

In the analysis given herein, which was made at the NACA Lewis labora- 
tory, the expression for eddy diffusivity given in reference 1 is 
modified in order to account for the effect of kinematic viscosity in 
reducing the turbulence in the region close to the wall. The effects 
of variable viscosity and of length-to-diameter ratio are also 
investigated. 


BASIC EQUATIONS 

‘For obtaining the velocity, temperature, and concentration distri- 
butions in a tube with turbulent flow, the differential equations for 
shear stress, heat transfer, and mass transfer can be written as fol- 
lows (symbols are defined in the appendix): 


du . du 

^ dy + pe d^F 

(1) 

dt dt 

dy " pgc P e h cLy 

(2) 

. dC dC 

dy h dy 

(3) 


where the values for e and are dependent on the amount and kind 

of turbulent mixing at a point. The eddy diffusivities for heat and 
mass transfer are equal inasmuch as both processes are governed by the 
same differential equation if aerodynamic heating is neglected (diffusion 
equation). On the other hand, the ratio e^/e = a must be determined 

experimentally or theoretically inasmuch as the equation of motion for 
a fluid contains terms which are not present in the diffusion equation. 
Equations (l) to (3) can be written in dimensionless form as 


t_ _ / _M_ + _P_ e \ du^ 

T o \^mo po pq/pq J d y + 


( 4 ) 
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JL- (— JL + A JjL. a —4 — \ (5) 

qo V*o Pr o pq c p,o Po/ p o/ciy 


_m_ _ / _X 1_ a e \ <3-C + 

mo “ \Xq Sc 0 Mq/pq/ dy + 


( 6 ) 


The variation of properties in these equations might be caused by 
either radial variation of temperature or of concentration of the 
diffusing substance. 


EXPRESSIONS FOR EDDY DIFFUSIVITY 


In order to make practical use of equations (4), (5), and (6), the 
eddy diffusivity e must be evaluated for each portion of the flow. 
After consideration of the various factors on which e might depend, 
the following functional relation is assumed: 


e 


e 



u du 

p' ay' 


d 2 u d 5 u 
dy 2 * dy 3 ' 


(7) 


(The quantities u and p must occur together inasmuch as they are 
the only quantities in eq. (7) containing mass as a dimension.) It is 
assumed in reference 8 that in the region at a distance from the wall 

e - e(du/dy, d 2 u/dy 2 ) (Khrman’s assumption), and in the region close to 
the wall 6 = e(u,y) = n 2 uy. In both of these expressions the possible 
effect of kinematic viscosity p/p is neglected. It appears from heat- 
and mass-transfer data at high Prandtl and Schmidt numbers, however, 
that the effect of p/p cannot be neglected in the region very close to 

the wall (y + < 5). When the previous expressions are used to compute 
heat and mass transfer at high Prandtl and Schmidt numbers, the coeffi- 
cients are considerably too high compared with the experimental values. 
It might be expected that, in the region very close to the wall where 
the turbulence level is low, the effect of kinematic viscosity would be 
important inasmuch as the viscous effects would be of the same order of 
magnitude as the inertia effects. Therefore, for the region close to 
the wall, e is written as 


6= e(u,y,p/p) 


( 8 ) 
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As in reference 8, the effect of the derivatives is neglected close to 
the wa.ll because the flow is very nearly laminar in that region. The 
first derivative approaches the value u/y and hence may be omitted 
since u and y already appear in the functional relation. The second 
derivative becomes very nearly zero as the wall is approached. 

From dimensional analysis, equation (8) becomes 

e = n^uyP . e'F^j (9, 


The function F[e'/(p/p)j should approach 1 as e or increases, 

because the effect of kinematic viscosity becomes negligible at high 
turbulence levels. Inasmuch as dimensional analysis cannot determine 
the form of the function F, additional assumptions must be made. 

The simplest assumption that might be made for F[e'/(p/p)] is 
that it equals e'/ ( m/p) * This assumption could be written in differen- 
tial form as 



However, equation (10) could not be expected to hold as F approaches 
1, because the change in F. for a given change in e'/(d/p) should 
approach 0 as F approaches 1 (F could never be greater than l). 

The simplest multiplicative factor which gives equation (10) this char- 
acteristic is (l - F). The assumption made for the variation of F, 
the adequacy of which will be checked by experiment, therefore 
becomes 


& = *[ e'/(n/p)J (l , f) 


(n) 


Separating variables and integrating from the wall to a point in the 
fluid give 


.F 


dF 


£ "/ (d/p) 


1 
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where the lower limits are taken as zero, because the turbulence goes 
to zero at the wall and the effect of kinematic viscosity consequently 
becomes extremely large at the wall (F *♦ 0). Integration of equation 
(12) results in 

' TmTpT 

F= 1 - e (13) 


This expression approaches 1 for large values of e 1 / (h/p) • From 
equations (9) and (13) there results 


e 




(14) 


Equation (14) gives e as a function of u and y for the region 
close to the wall. ' The constant n is to he determined experimentally. 

For the region at a distance from the wall, F is usually close to 
1, because the effect of kinematic viscosity is small. The Karman 
expression for which neglects the variation of F, can usually be 
used in that region, or 


e - x 2 (du/dy) 5 
( d 2 u/dy 2 ) 2 

If it is desired to take the variation of F into account in the region 
at a distance from the wall, equation (15) becomes 


x 2 (du/dyf/ (d 2 u/dy 2 ) 2 
IVP 



e = x 2 


(du/dy) 5 

(d 2 u/dy 2 ) 2 


1 - e 



No attempt is made in the present analysis to specify the mechanism 
by which the kinematic viscosity reduces the eddy diffusivity in the 
region close to the wall, because the exact mechanism is unknown. Pos- 
sible mechanisms are these: First, the kinematic viscosity might act as 

a dancing factor to reduce the turbulence level close to the wall. 
Second, it might help to orient the eddies close to the wall by damp- 
ing. out those moving perpendicular to the wall and thus reduce tne 
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effective turbulent transfer. Third, it might act to produce a partial 
turbulence; that is, it might cause the flow at a point to be laminar a 
fraction of the time. The actual effect of kinematic viscosity might 
be due to any, or all, of these mechanisms. 


ANALYSIS FOR CONSTANT FLUID PROPERTIES 

In order to solve equations (4), (5), and (6), the following 
assumptions are made in addition to the assumptions concerning the 
expressions for eddy diffusivity (eqs. (14), (15), and (16)): 

(1) The eddy diffusivities for momentum e and heat or mass trans- 
fer Cfc are equal, or a = 1. Previous analyses for flow of gases in 
tubes based on this assumption yielded heat-transfer coefficients that 
agree with experiment (ref. l). At low Peclet numbers (Pe = RePr), a 
appears to be a function of Peclet number (ref. 9) but is approximately 

1 at high Peclet numbers. In general, the Peclet numbers are high at high 
Prandtl numbers for turbulent flow. 

(2) The variations across the tube or boundary layers of shear 
stress t , the heat transfer per unit area q, and the mass transfer per 
unit area m have a negligible effect on the velocity, temperature, and 
concentration distributions. It is shown in figure 12 of reference 1 
that the assumption of a linear variation of shear stress and heat trans- 
fer across the boundary layers gives very nearly the same velocity and 
temperature profiles for gases as those obtained by assuming uniform 
shear stress and heat transfer across the boundary layers for values 

of 5^ + between 500 and 5000. (The boundary layers fill the tube for 
fully developed flow or fully developed heat or mass transfer.) For 
small values of 5h + > such as occur very near the entrance, the effect 
of variable heat transfer (or mass transfer) is checked in figure 9 of 
reference 10 and found to be negligible. Although these checks were 
made for gases (Prandtl or Schmidt numbers close to l), the effect of 
variation of heat or mass transfer per unit area at high Prandtl or 
Schmidt numbers would be even less because the temperature or concen- 
tration profile becomes flatter as the Prandtl or Schmidt number is 
increased. 

(3) The molecular shear-stress, heat -transfer, and mass-transfer 
terms in the equations can be neglected in the region at a distance from 
the wall (ref. 1, fig. 14). 

(4) In the case of mass transfer, the concentration of the dif- 
fusing substance is small enough that the mass transfer does not appreci- 
ably change the velocity. This condition is generally obtained in the 
case of evaporation from a wetted wall or in the solution of the wall 
material in a liquid. 
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Velocity, temperature, and, concentration distributions . - The ex- 
pression for e close to the wall (eq. (14)) can be written in dimen- 
sionless form for constant properties as 



n 2 u +y+ 



e 


n 2 u + y +> 


fl7) 


Equations (4), (5), and (6) can be written in integral form for the 
region close to the wall with the preceding assumptions and constant 
fluid properties as 



(18) 



(19) 



( 20 ) 


For the region at a distance from the wall, equation (15) is substituted 
in equations (4), (5), and (6). By use of assumptions (2) and (3), 
equation (4) becomes, for constant fluid properties. 


< . i lo ge ( ) ( 2 1) 

which is the well-known logarithmic equation. Division of equation (4) 
by equation (5) or (6) gives, with assumptions (l), (2), and (3), 
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u + - = t + - t-j* = C + - C x + (22) 

where the equations are integrated from y^ + , the lowest value of y + 
for which the equations for flow at a distance from the wall apply, to 


The variation of u + with y + for fully developed adiabatic tur- 
bulent flow from references 8 and 11 is shown on semilogarithmic co- 
ordinates in figure 1. The curves corresponding to equations (18) and 
(21) are also shown in the figure. Equation (18) was solved by numer- 
ical iteration inasmuch as u + occurs on both sides of the equation. 

The values of the constants in the equations are n = 0.124, x = 0.36, 

and y^ + = 26, as determined from the experimental data. These values 

will apply also to the temperature distributions. The value of n in 
equation (l6) differs from that given in reference 8 because equation 

(18) includes the effect of kinematic viscosity on e* Figure 1 indi- 
cates good agreement of equation (18) with the data for y+ < 26 and of 
equation (2l) with the data for y + > 26. From the velocity-distribution 

data it is difficult to tell whether e = n^uy from reference 8 or 
equation (14) for e should be used in the region close to the wall, 
inasmuch as both expressions give results which agree closely with the 
data. The temperature or concentration profiles at high Prandtl or 
Schmidt numbers are, however, much more sensitive to the values of e 
very close to the wall because of the very large temperature or concen- 
tration gradients in that region, as can be seen in figure 2. 

Generalized temperature or concentration distributions calculated 
from equations (19), (20), and (22) are presented in figure 2 on log- 
log coordinates. Each curve represents either t + or C + at a given 
Prandtl or Schmidt number, as can be seen by comparison of equations 

(19) and (20). The curves indicate that the temperature or concentration 
distributions become flatter over most of the tube radius as the Prandtl 
or Schmidt number increases. From equations (19) and (20), dt + /dy + = Pr 
and dC + /dy + = Sc at or very near the wall so that the slopes of the 
curves at the wall increase with Prandtl or Schmidt number. The slopes 
of the curves in figure 2 near the wall appear equal because the curves 
are plotted on log-log coordinates (d(log t+)/d(log y+) = 1 at the wall). 
Included for comparison is the temperature distribution for a Prandtl 
number of 300 calculated by the method in reference 1, which neglects 

the effect of kinematic viscosity (e = n^uy). 

The sensitivity of the temperature or concentration distribution at 
high Prandtl or Schmidt numbers to various assumptions for the turbulent 
transfer in the region close to the wall compared with that of the 
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velocity distribution indicates that the region very close to the wall 
could he studied advantageously hy measuring temperatures or concen- 
trations at high Prandtl or Schmidt numbers rather than by measuring 
velocities in that region. Some work along these lines has been reported 
in reference 5, in which concentration profiles at high Schmidt numbers 
were measured with an interferometer. No evidence of a purely laminar 
layer (linear concentration profile) was found for values of y + as 
low as 1. This result is in agreement with the assumption in the pre- 
sent analysis, where the turbulence is assumed to be zero only at the 
wall. 


Relations among Nusselt , Reynolds , and Prandtl or Schmidt numbers 
for constant properties . It can be shown from the definitions of the 
quantities involved that the Nusselt numbers for heat and mass transfer 
and the Reynolds number are given by 


Nu = 


2rp + Pr 

*b + 


(23) 


Nu' 


2r Q + Sc 
“^ 5 ~ 


(24) 


Re = 2u b + r 0 + 


(25) 


vhere 



(26) 



( 27 ) 
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and 



(28) 


The Nusselt numbers in these equations are based on the difference 
between wall and bulk temperatures or concentrations. The relation 
among Nusselt, Reynolds, and Prandtl numbers can be obtained from these 
equations and the generalized distributions given in figures 1 and 2 
for various values of the parameter r Q + . 

Predicted Nusselt numbers for fully developed heat or mass transfer 
are plotted against Reynolds number for various values of Prandtl or 
Schmidt n umb er in figure 3. The curve for a Prandtl number of 0.73 
agrees very closely with that given in reference 1, which was obtained 

by using e = n^uy in the region close to the wall. The curves for 
supercritical water at higher Prandtl numbers (Pr » 6) given in refer- 
ences 12 (fig. 7) and 13 (fig. 18), however, are higher than those in 
figure 3; therefore the curves in these references should be replaced 
by figure 3, although the values of the reference temperatures for 
evaluating the fluid properties should not be affected. 

Examination of the curves in figure 3 indicates that the slopes of 
the various curves are approximately equal on a log-log plot . (The 
slopes would be more nearly equal if the Peclet number effect from 
ref. 9 were included.) This result justifies the usual practice in 
heat-transfer investigations of writing Nu = f(Re,Pr) as f(Re)x f(Pr) 
(usually as Re a Pr b ) . The same result does not hold for very low 
Prandtl numbers where the slopes change considerably. 

Comparison of analysis and experiment for fully developed heat and 
mass transfer. - A comparison between predicted and experimental results 
for fully developed heat and mass transfer is presented in figure 4, in 
which Stanton’ number is plotted against Prandtl or Schmidt number for 
various Reynolds numbers. The predicted Stanton numbers were obtained 
from figure 3 and the relation St = Nu/RePr. The symbols represent 
mean lines through data for heat transfer in gases (ref. l) and in 
liquids (refs. 14 to 19) and mass transfer by evaporation from wetted 
walls (refs. 20 to 22), by solution of the wall material in a liquid (refs 
23 and 24), and by diffusion-controlled electrodes (ref. 25). The pre- 
dicted and measured values are in good agreement over the entire range 
of Prandtl and Schmidt numbers shown (0.5 to 3000). The agreement for a 
Reynolds number of 10,000 in the low Prandtl or Schmidt number range 
would be improved by applying the Peclet number correction from 
reference 9. 
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Simplified equation for fully developed case for very high Pratidtl 
or Schmidt numbers. - At very high Prandtl numbers the essential tem- 
perature changes take place in the region very close to the wall where 

u + is very nearly equal to y + . Setting u + = y + in equation (19), 
expanding the exponential function in a series, and retaining only the 
first two terms of the series result in 



Integration of equation (29) and evaluation of the result for 
give 


(29) 




it 

2/^2n 



(30) 


The value of t + at y + = «. is essentially % + , because t + is very 

nearly constant except in the region very close to the wall at high 
Prandtl numbers. From equations (23) and (30), 


Nu = nrQ + PrV 4 (3l) 

Jt 

or the Stanton number, in terms of the Prandtl number and friction 
factor , is 


St 


2n 'V f~ 
a p r 3/ 4 


where 



Similarly, for mass transfer, 


(32) 


St 


I 


2n VF 
* Sc 3 / 4 


(33) 


where n has the value 0.124 as determined in figure 1. The relation 
between Re and f is given in figure 8 of reference 1 or the curve 
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in figure 10 in the present report for adiabatic flow (0=0). This 
relation is, of course, independent of Prandtl number for constant 
properties. Equation (32) or (33) is indicated by the dotted line in 
figure 4 and is seen to be in good agreement with the predicted line 
obtained previously for Pr>200. 

Comparison of various analyses . - A comparison of various analyses 
is given in figure 5. It can be seen that all the analyses- more or less 
converge at the lower Prandtl or Schmidt numbers. At the high Prandtl 
or Schmidt numbers, the present analysis and the analyses from refer- 
ences 4 and 5 are in fair agreement, whereas those from references 2 
and 3 diverge considerably. The present analysis and the analysis from 
reference 5 represent the experimental data about equally well. 

Heat or mass transfer in entrance region for uniform wall heat or 
mass flux, uniform initial temperature or concentration distribution, 
fully developed velocity distribution, and constant properties . - For 
calculating heat or mass transfer in the entrance region, it is assumed, 
as in reference 10, that the effects of heat or mass transfer are con- 
fined to fluid layers close to the surface (boundary layers for heat or 
mass transfer). The temperature or concentration distributions outside 
the boundary layers are assumed uniform, and the temperature or concen- 
tration is constant along the length of the tube for the region outside 
the boundary layer. Inside the boundary layer the temperature or con- 
centration distribution is obtained from figure 2. Integral heat- or 
mass -transfer equations are used for calculating the thickness of the 
boundary layers for various distances from the entrance. It is shown in 
reference 10 that for heat transfer the relation between 8^ and x/D 

for constant properties is given by 


X 

D 


+ ,2 


2(0 


(t 5 + - t + )u + (r Q + - y + )dy + 


(34) 


'0 


A similar equation can be written for mass transfer if the concentration 
of the diffusing material is small (assumption 4): 


V 


X 

D 


+ \2 


2 ( r O ) 


(C 5 + - C + )u + (r Q + - y + )dy + 


'0 


(35) 
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In equations (34) and (35), X represents the axial distance from the 
point at which heat or mass transfer begins. The dimensionless boundary- 
layer thickness &k + is the same in both equations when the Prandtl 
number equals the Schmidt number inasmuch as C* = t + for a given value 
of y*" and Prandtl or Schmidt number. 

Values of local Nusselt and Reynolds numbers can be obtained from 
equations (23), (24), and (25), as for fully developed flow, with the 
exception that the expression for tt> + in equation (26) is replaced by 



/ \ 

A similar expression can be obtained for Cv, ^ The integral in the 
v : • ' j , w 

numerator is broken into two parts, because v is constant and eq ua l 

to t$ + outside the thermal boundary layer. Inside the thermal bound- 
ary layer the relation between t + and y* is obtained from figure 2. 
Values of u + are obtained from figure 1 for y + from 0 to Tq + 

inasmuch as a fully developed velocity distribution is assumed. The 
relation between Nusselt number and x/D for various values of Reynolds 
number is obtained by assuming values of the parameters Tq + and 5^+ 
and by calculating the various quantities from equations (23), (25), 

(35), and (36). 

The variation of local Nusselt number for heat or mass transfer 
divided by the fully developed Nusselt number with x/D and Reynolds 
number at Prandtl or Schmidt numbers between 1 and 3000 is given in 
figure 6. At the higher Reynolds numbers the values of Nu/Nu^ for a 

given x/D decrease with increasing Prandtl number; that is, the effect 
of x/D becomes small at large values of Prandtl number. At low Rey- 
nolds numbers the variation is more complex: Values of Nu/Nu^ first 

decrease and then increase slightly as Prandtl number increases. In 
either case it is evident that in the entrance region the fully developed 
Nusselt numbers should be multiplied by a factor which is a function 
of X/D, Reynolds number, and Prandtl number rather than of x/D alone. 
That is, a simiple factor such as (x/d)“ c , which is often used for Prandtl 
numbers of approximately 1, is inadequate for high Prandtl numbers. 

The effect of Reynolds number on Nusselt number in the entrance region 
increases; that is, the separation of the curves for various Reynolds 
numbers increases with Prandtl number. The same conclusions apply to 
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the average Nusselt numbers plotted in figure 7 except that the changes 
with x/D near the entrance are more gradual and the separation of the 
curves with Reynolds number is greater than for the local values. 

The average Nusselt numbers were calculated from 



This equation is consistent with the definition 


(*0 " t b) a v = 



(t 0 - t^dX 

~X 


which is the usual way of defining the average difference between wall 
and bulk temperatures for uniform heat flux. 

It is of interest to note that the Nusselt numbers for turbulent 
flow in figures 6 and 7, in general, display trends opposite to those 
for Nusselt numbers for laminar flow with increasing values of Prandtl 
number. In the case of laminar flow, the value of Nu/Nu^ at a given 

X/D and Reynolds number increases considerably with Prandtl number 
(ref. 26, fig. 3), because the heat diffuses through the fluid more 
slowly at the higher Prandtl numbers (the thermal diffusivity is lower) 
so that the thermal boundary layer is thinner and Nu/Nu^ for a given 
X/D near the entrance and a given Reynolds number is consequently 
higher than for the lower Prandtl numbers. The same phenomenon also 
tends to increase the effect of x/D for turbulent heat transfer. In 
the case of turbulent heat transfer, however, the shape of the tem- 
perature profile in the thermal boundary layer changes considerably with 
Prandtl number and becomes very flat at high Prandtl numbers (fig. 2). 
This means that the temperature profiles for fully developed flow do 
not differ greatly from those near the entrance (both are flat) although 
the boundary-layer thickness for the two cases differs considerably. 

The Nusselt numbers in the entrance region for turbulent heat transfer 
at high Prandtl numbers therefore tend to quickly approach the fully 
developed values. 
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FULLf DEVELOPED HEAT TRANSFER AND FRICTION WITH VARIABLE PROPERTIES 

Turbulent heat transfer to liquids with variable viscosity . - In 
the case of heat transfer to liquids, the variation of the viscosity 
with temperature is considerably greater than the variation of the 
other properties. A good approximation to the actual heat transfer in 
liquids can therefore be obtained by considering only the viscosity to 
vary with temperature. Under that assumption and assumptions 1 and 2 
from the section Analysis for Constant Fluid Properties, equations (4) 
and (5) become 


and 



du + 

dy+ 

dt + 

dy + 


(37) 


(38) 


For variable viscosity, the expression for e close to the wall (eq. 
(14)) can be written in dimensionless form as 


^o/P 0 


n 2 u + y + 



n 2 u + y + 

p/pq 


(39) 


For liquids, including water, oil, ethylene glycol, and sodium hydroxide, 
h/Fq can usually be represented by (t/tQ), if the liquid is noa; too near 
the freezing point. The exponent d varies from -1 to -4, and the temper- 
atures are measured in °F. This differs from the case for gases where the 
temperatures were measured in °R (ref. l). All the results up to this 
point are independent of whether the temperatures are in °R or °F. From 
the definitions of t + and p, 

li/liQ = (1 - Pt + ) d (40) 

Substituting equations (39) and (40) into equations (37) and (38) and 
writing the result in integral form yield 
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and 



Equations (4l) and (42) can be solved simultaneously by iteration, that 
is, assumed values for u+, y + , and t + are substituted into the right 
side of the equations and new values of u + and t + are calculated by 
numerical integration. These new values are then substituted into the 
right side of the equations and the process is repeated until the values 

of u + and t + corresponding to each value of y + do not change 
appreciably. Equations (4l) and (42) give the relations between u + , 
t , and y + for various values of the heat -transfer parameter 3 and 
of Prg for the region close to the wall. 

In the region at a distance from the wall, the terms in the equa- 
tions containing variable viscosity are neglected so that u + and t + 
are given by equations (21) and (22). 

Generalized temperature and velocity distributions for a Prandtl 
number at the wall of 10 and a d of -4 are shown for various values 
of the heat -transfer parameter 3 in figure 8. The value y + at the 
intersection of the curves for flow close to and at a distance from the 

wall is taken at y^* = 26, as in the case of constant viscosity. The 
effect of various assumptions for the variation of y-j* is investigated 
in reference 13 (fig. 13), in which it is concluded that this assumption 
should give accurate results. Positive values of 3 correspond to heat 
addition to the liquid, negative values to heat extraction. The values 

of t + at a given value of y+ increase with increasing 3> whereas 
the values of u + decrease. These opposite trends can be explained 
by examining the terms in the denominators of equations (4l) and (42). 
The first term in the denominator of equation (41 ) (molecular) causes 
u + to decrease with increasing values of 3 (d < 0), whereas the 
second term (turbulent) causes u + to increase. The effect of the 
second term is somewhat smaller so that the net effect is a decrease, as 
shown in figure 8(b). In equation (42), however, the first term is 
independent of 3 and therefore the second term causes an increase in 
t + with increasing values of 3. The opposite trends for heat transfer 
and friction in most of the succeeding curves can be attributed to the 
same cause. 
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The effect of the factor in parenthesis in equation (42), which is 
equal to F in equation (13) , increases with increase in Prandtl number 
because of the steep temperature gradients at high Prandtl numbers in 
the region very close to the wall where F differs considerably from 1. 
The separation of the t + curves for various values of p (not shown) 
therefore increases as Prandtl number increases, the separation for a 
Prandtl number of 1 being very small. 

For obtaining the relation between Nusselt number, Reynolds number, 
Prandtl number, and friction factor for variable viscosity, equations 
(23), (25), and (26) apply for variable properties as well as for con- 
stant properties if the viscosities in the Reynolds and Prandtl numbers 
are evaluated at the wall temperature . The friction factor f can be 
calculated from 


f = 2/(u b +) 2 (43) 

Nusselt numbers and friction factors are plotted against Reynolds 
numbers in figures 9 and 10, respectively, with the viscosity in the 
Reynolds and Prandtl numbers evaluated at the wall temperature. The 
curves for cooling and heating are for values of p b /pQ on the order of 

0.5 and 2, respectively. As in the case of temperature and velocity 
distributions, the trends with increasing values of p for the Nusselt 
numbers and friction factors are opposite. Also, as in the case of the 
temperature distributions, the separation of the Nusselt number curves 
with p increases as the Prandtl number increases, the separation being 
very small for a Prandtl number of 1. Thus, for a Prandtl number of 1 
the reference temperature for evaluating the viscosity in order to 
eliminate the effects of variable viscosity is close to the wall tem- 
perature; the departure of the reference temperature from the wall tem- 
perature increases with Prandtl number. This is to be expected because, 
as mentioned previously in this section, the separation of the t + 
against y+ curves with p increases with Prandtl number. In those 
curves all the properties are evaluated at the wall temperature. 

The reference temperature for heat transfer in liquid with variable 
viscosity for a Prandtl number of 1 differs from that obtained for 
gases in reference 1, in which all the properties except the specific 
heat were considered variable. In reference 1, the reference temperature 
for gases was found to be close to the average of the wall and bulk 
temperatures rather than close to the wall temperature. In assigning a 
reference temperature, it is therefore important to consider what 
properties are variable. 

Values of x for calculating the reference temperature t x , where 
t x = x(t 0 - t b ) + tfc, are given, for heat transfer and friction, in 
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figure 11. The curves for 3=0 in figures 3 and 10 can "be used for 
variable viscosity if the viscosities in the Reynolds and Prandtl num- 
bers are evaluated at the reference temperatures given in figure 11 
(Re x and Pr x ) . The values of x were computed for values of d 
(p/|Jq = (t/tQ)^) of -1 and -4 and for values of M-^/^o a ^ ou ^ 0*5 

and 2. The value of d had little effect on the curves, but different 
curves are obtained for heating and cooling of the liquid. In the case 
of heat transfer, the reference temperature does not depart greatly from 
that in the widely used Colburn equation (ref. 27), wherein the viscosity 


is evaluated at t, 


0.5 


except at the lower Prandtl numbers. The values 


of x for friction are lower than those for heat transfer. Deviations 
from the curves in figure 11 might occur for very high viscosity ratios 
or for cases in which the viscosity variation with temperature could 
not be represented by a simple power function. 


Turbulent heat transfer to gases with variable properties using 
present method of analysis . - If the present method of analysis is to 
be considered more general than the analysis in reference 1, it must 
be applicable to gases with variable properties as well as to liquids. 
When assumptions (l) and (2) in the section Analysis for Constant Fluid 
Properties are used, equations (4) and (5) become, for gases. 


and 


+ _P_ e \ du + 

^0 p o p o/ p O / dy + 


k 1 t p e \ dt + 

*0 ^0 + P 0 V p 0 ) dy + 


(44) 


(45) 


where Cp is assumed constant "because its variation with temperature is 
slight compared with the variations of viscosity, thermal conductivity, 
and density. As in reference 1, it is assumed that k/k^ = 
and p/pQ = tQ/t, where d = 0.68 and the temperatures are measured in 
degrees Rankine. By substituting the expression for e close to the 
wall (eq. (14)) and t/t q = 1 - pt + , equations (44) and (45) can be 
written in integral form for the region close to the wall as 



n 2 u + y + \ 

(i - ut + ) d+1 j 


( 46 ) 
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dy + 



n 2 u + y + 


(1 - (3t + ) 


d+1 


(47) 


Equations (46) and {47 ) can be solved simultaneously by iteration as 
were equations (4l) and (42). 

In the region at a distance from the wall, the molecular shear- 
stress and heat-transfer terms in equations (44) and (45) are neglected 
and the expression for e given in equation (15) is used. The inte- 
gration is carried out in reference 1, in which it is found that 



where y^ + = 26. 

Generalized temperature and velocity distributions as calculated by 
the present analysis for a Prandtl number of 0.73 are plotted in figure 
12. The distributions calculated in reference 1 are also included for 
comparison. The agreement between the two methods of analysis is 
satisfactory, so that the Nusselt numbers and the reference temperatures 
obtained from the present analysis should also agree with those given 
in reference 1. 

Laminar heat transfer to liquids with variable viscosity and 
uniform heat flux . - Laminar heat transfer to liquid metals with variable 
viscosity is investigated in reference 28. The results of that inves- 
tigation should be applicable also to liquids with high Prandtl numbers 
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if the value of d is on the order of -1.6. The temperatures in the 
report are to be measured in °F (m/pq = (t/tQ)^). It is found (ref. 
28) that the results for fully developed laminar flow with constant 
heat flux can be represented closely by 

- (8)(M° “ 


SUMMARY OF RESULTS 

The following results were obtained from the analytical investi- 
gation of heat and mass transfer in smooth tubes at high Prandtl and 
Schmidt numbers: 

1. By modifying the expression for eddy diffusivity from a previous 
analysis to account for the effect of kinematic viscosity in the region 
close to the wall, good agreement was obtained between predicted and 
experimental results for heat and mass transfer at Prandtl and Schmidt 
numbers between 0,5 and 3000. A simplified equation was obtained for 
very high Prandtl or Schmidt numbers. 

2. The analysis indicated that, except at low Reynolds numbers, 
the entrance effect (local Nusselt number divided by fully developed 
Nusselt number) for heat or mass transfer decreases as Prandtl or 
Schmidt number increases. 

3. The analysis indicated that the effects of variable viscosity 
on turbulent heat transfer and friction in liquids can be nearly 
eliminated in ordinary cases by evaluating the viscosities in the Rey- 
nolds and Prandtl numbers at reference temperatures which are functions 
of the Prandtl number. For the laminar case the results for liquid 
metals with variable viscosity given in a previous analysis should be 
applicable to liquids at high Prandtl numbers. 

4. When the present method of analysis was applied to gases with 
variable properties, essentially the same results were obtained as are 
reported from analysis and experiments in a previous report. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, February 17, 1954 
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APPENDIX - SYMBOLS 

The following symbols are used in this report: 

C concentration of diffusing substance, (lb) (sec 2 )/ft 4 

C b bulk concentration of diffusing substance, (lb) (sec 2 )/ft 4 

Cg concentration of diffusing substance at y = 5, 

(lb) (sec 2 )/ft 4 

Cq concentration of diffusing substance at wall, 

(lb) (sec 2 )/ft 4 

c exponent 

c p specific heat of fluid at constant pressure, Btu/(lb)(°F) 

c D q specific heat of fluid at constant pressure at wall, 

' Btu/(lb)(°F) 

D inside diameter of tube, ft 

d exponent, value of which depends on variation of viscosity 

of fluid with temperature 

F function of e'/(p/p) 

g acceleration due to gravity, 32.2 ft/sec^ 

h local heat-transfer coefficient, Qo/(to - t b ), 

Btu/ (sec) (sq ft) (°F) 

h' local mass-transfer coefficient, n^/(C 0 - C b ), ft/sec 

h av average heat-transfer coefficient, %/(t 0 - t-fcJav 

h av ' average mass-transfer coefficient, n^/(t 0 - t b ) 

k thermal conductivity of fluid, Btu/ (sec) (sq ft)(°F/ft) 

kQ thermal conductivity of fluid evaluated at tn, 

Btu/(sec)(sq ft)(°F/ft) 

m rate of mass transfer toward tube center per unit area, 

(lb) (sec)/cu ft 
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n 


% 

r 0 

t 

*b 


0 

U 

% 

x 


■ tb) 


av 


y 

&h 

e 

e • 


rate of mass transfer toward tube center per unit area 
at wall, (lb)(sec)/cu ft 

constant 

rate of heat transfer toward tube center per unit area, 
Btu/ (sec) (sq ft) 

rate of heat transfer at inside wall toward tube center 
per unit area, Btu/ (sec) (sq ft) 

inside tube radius, ft 

temperature, °F for liquids or °B for gases 

bulk static temperature of fluid at cross section of tube, 
°F 

reference temperature for local Reynolds and Prandtl 
numbers, x(t Q - t t ) + t^, °p 

temperature of fluid outside thermal boundary layer, °F 

wall temperature, °F for liquids or °R for gases 

average difference between wall and bulk temperature, °F 

time-average velocity parallel to axis of tube, ft/sec 

bulk velocity at cross section of tube, ft/sec 

distance from point at which heat or mass transfer 
begins, ft 

number used in evaluating arbitrary temperature in tube 
t x. 

distance from wall, ft 

thermal or diffusion boundary -layer thickness, ft 

coefficient of eddy diffusivity for momentum, sq ft/sec 

expression for eddy diffusivity which neglects effect of 
l-t/ p, sq ft/sec 


coefficient of eddy diffusivity for heat or mass, sq ft/sec 
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X 

X 0 

h 

^b 

M- 

x 

p 0 

P 

P 0 

x 

T 0 


Karman constant 

molecular diffusivity, sq ft/sec 

molecular diffusivity at wall, sq ft/sec 

absolute viscosity of fluid, (lb)(sec)/sq ft 

absolute viscosity of fluid evaluated at t , 
(lb)(sec)/sq ft 

absolute viscosity of fluid evaluated at t , 
(lb)(sec)/sq ft X 

absolute viscosity of fluid evaluated at tQ, 
(lb)(sec)/sq ft 

mass density of fluid, (lb) (sec^/ft 4 

mass density of fluid evaluated at tQ, (lb) (sec^)/ft^ 

shear stress in fluid, lb/sq ft 

shear stress in fluid at wall, lb/sq ft 


Dimensionless groups: 


Cb 


C 5 + 


C 

°1 

f 

Nu 

Nu’ 


Nu 


av 


Nu. 


concentration parameter. 


( c o ~ c )V t q/pq 

“b 

bulk concentration parameter, (Cq - c b )V T o/ p 0 

“b 


value of C + at 5 V 


av 


+ , ( C 0 " C 5 ^ / T o/ p 0 

n 7 

“b 

value of C + at y^ + 
friction factor, 2 Tq/pu^ 

Nusselt number for heat transfer, hD/k 
Nusselt number for mass transfer, h'D/\ 
average Nusselt number for heat transfer, h &v D/k 
average Nusselt number for mass transfer, h av ' D A 
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Nu d 

fully developed Wusselt number for heat transfer 


Nu d ’ 

fully developed Nusselt number for mass transfer 


Pe 

Peclet number, pu^Dcpg/k 


Pr 

Prandtl number, c p gp/k 


Pr 0 

Prandtl number with properties evaluated at tQ 


Re 

Reynolds number, pu^D/p 


Re x 

Reynolds number with viscosity evaluated at t x . 

P u b D /^ x 

Re 0 

Reynolds number with viscosity evaluated at tn> 

P u b D /Mo 

r 0 + 

MVpo 

tube radius parameter, Vq/p r o 

Sc 

Schmidt number, p/(p\) 


Sc 0 

Schmidt number at wall, Po/po\> 


St 

Stanton number for heat transfer, h/pgu^c^ 


St’ 

Stanton number for mass transfer, h'/ u d 
temoerature Darameter. — 2 — ^ ^ ® 


t + 


V 

i 

bulk-temperature parameter, i | 1 - j 

i/ M 1 0> 



p V 1 ' *0 / 


h + 

value of t + at y^ + 


u + 

velocity parameter, wj t q/Pq 


u b + 

bulk -velocity parameter, u ^^q/Pq 


u l + 

value of u + at y-^ + 

A f T 0 / Pn 

wall distance parameter, — — , y 


y + 
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*1 


value of y + at intersection of curves for flow close to 
' wall and at a distance from wall 


a 

3 


ratio of eddy diffusivities, e h /e 

heat -transfer parameter, q 0 V T o/ P o/ ( c p^0 t 0 ) 

dimensionless thermal or diffusion boundary -layer thick- 

W p 0 

ness, -7 — — 8 h 

^0/P 0 h 


REFERENCES 

1. Deissler, R. G., and Eian, C. S.: Analytical and Experimental 

Investigation of Fully Developed Turbulent Flow of Air in a Smooth 
Tube with Heat Transfer with Variable Fluid Properties. NACA 
TN 2629, 1952. 

2. von Karman, Th.: The Analogy Between Fluid Friction and Heat Transfer. 

Trans. A.S.M.E., vol. 61, no. 8, Nov. 1939, pp. 705-710. 

3. Summerfield, Martin: Recent Developments in Convective Heat Transfer 

with Special Reference to High-Temperature Combustion Chambers. 

Heat Transfer Symposium, Eng. Res. Inst., Univ. Mich., 1953, 
pp. 151-171. 

4. Murphree, E. V.: Relation Between Heat Transfer and Fluid Friction. 

Ind. and Eng. Chem. , vol. 24, no. 7, July 1932, pp. 726-736. 

5. Lin, C. S., Moulton, R. W. , and Putnam, G. L.: Mass Transfer Between 

Solid Wall and Fluid Streams. Ind. and Eng. Chem., vol. 45, no. 3, 
Mar. 1953, pp. 636-640. 

6. Sherwood, Thomas K. : Heat Transfer, Mass Transfer, and Fluid 

Friction - Relationships in Turbulent Flow. Ind. and Eng. Chem., 
vol. 42, no. 10, Oct. 1950, pp. 2077-2084. 

7. Martinelli, R. C.: Heat Transfer to Molten Metals. Trans. A.S.M.E., 

vol. 69, no. 8, Nov. 1947, pp. 947-959. 

8. Deissler, Robert G.: Analytical and Experimental Investigation of 

Adiabatic Turbulent Flow in Smooth Tubes. NACA TN 2138, 1950. 

9. Deissler, Robert G.: Analysis of Fully Developed Turbulent Heat 

Transfer at Low Peclet Numbers in Smooth Tubes with Application to 
Liquid Metals. NACA RM E52F05, 1952, 



26 


NACA TN 3145 


10. Deissler, Robert G.: Analysis of Turbulent Heat Transfer and Flow 

in the Entrance Regions of Smooth Passages. NACA TN 3016, 1953. 

11. Laufer, John: The Structure of Turbulence in Fully Developed Pipe 

Flow. NACA TN 2954, 1953. 

12. Deissler, Robert G., and Taylor, Maynard F.: Analysis of Heat 

Transfer and Fluid Friction for Fully Developed Turbulent Flow of 
Supercritical Water with Variable Fluid Properties in a Smooth 
Tube. NACA RM E53B17, 1953. 

13. Deissler, R. G.: Heat Transfer and Fluid Friction for Fully Devel- 

oped Turbulent Flow of Air and Supercritical Water with Variable 
Fluid Properties. Trans. A.S.M.E., vol. 76, no. 1, Jan. 1954. 

14. Bernardo, Everett, and Eian, Carroll S.: Heat Transfer Tests of 

Aqueous Ethylene Glycol Solutions in an Electrically Heated Tube. 
NACA WR E-136, 1945. (Supersedes NACA ARR E5F07.) 

15. Kaufman, Samuel J., and Isely, Francis D.: Preliminary Investigation 

of Heat Transfer to Water Flowing in an Electrically Heated 
Inconel Tube. NACA RM E50G31, 1950. 

16. Eagle, A. E., and Ferguson, R. M. : On the Coefficient of Heat 

Transfer from the Internal Surface of Tube Walls. Proc. Roy. Soc. 
(London), ser. A, vol. 127, 1930, pp. 540-566. 

17. Kreith, F., and Summerfield, M. : Pressure Drop and Convective Heat 

Transfer with Surface Boiling at High Heat Flux; Data for Aniline 
and n -Butyl Alcohol. Trans. A.S.M.E., vol. 72, no. 6, Aug. 1950, 
pp. 869-878; discussion, pp. 878-879. 

18. Grele, Milton D., and Gedeon, Louis: Forced-Convection Heat-Transfer 

Characteristics of Molten Sodium Hydroxide. NACA RM E52L09, 1953. 

19. Hoffman, H. W.: Turbulent Forced Convection Heat Transfer in Cir- 

cular Tubes Containing Molten Sodium Hydroxide. 0RNL-1370, 

Reactor Exp. Eng. Div., Oak Ridge Nat. Lab., Oak Ridge (Tenn.), 

Oct. 3, 1952. (Contract No. W-7405, eng. 26.) 

20. Barnet, Walter X., and Kobe, Kenneth A.: Heat and Vapor Transfer 

in a Wetted-Wall Tower. Ind. and Eng. Chem. , vol. 33, no. 4, 

Apr. 1941, pp. 436-442. 

21. Chilton, T. H., 'and Colburn, A. P.: Mass Transfer (Absorption) 

Coefficients. Ind. and Eng. Chem., vol. 26, no. 11, Nov. 1934, 
pp. 1183-1187. 



NACA TN 3145 


27 


22. Jackson, M. L., and Ceaglske, N. H. : Distillation, Vaporization, 

and Gas Absorption in a Wetted-Wall Column. Ind. and Eng. Chem. , 
vol. 42, no. 6, June 1950, pp. 1188-1198. 

23. Bonilla, Charles F.: Mass Transfer in Liquid Metal and Fused Salt 

Systems. NY0-3086, First Quarterly Prog. Rep., U. S. Atomic 
Energy Comm., Tech. Information Service, Oak Ridge (Tenn.), 

Sept. 1, 1951. (Contract No. AT(30-l)-1100. ) 

24. Linton, W. H., Jr., and Sherwood, T. K. : Mass Transfer from Solid 

Shapes to Water in Streamline and Turbulent Flow, Chem. Eng. Prog 
vol. 46, no. 5, May 1950, pp. 258-264. 

25. Lin, C. S. , Denton, E. B., Gaskill, H. S., and Putnam, G. L. : 

Diffusion-Controlled Electrode Reactions. Ind. and Eng. Chem., 
vol. 43, no. 9, Sept. 1951, pp. 2136-2143. 

26. Poppendiek, H. F., and Palmer, L. D.: Forced Convection Heat Trans 

fer in Thermal Entrance Regions, Pt. II. ORNL 914, Metallurgy 
and Ceramics, Reactor Exp. Eng. Div., Oak Ridge Nat. Lab., Oak 
Ridge (Tenn.), May 26, 1952. (Contract No. W-7405, eng. 26.) 

27. McAdams, William H.: Heat Transmission. Second ed., McGraw-Hill 

Book Co., Inc., 1942, p. 168. 

28. Deissler, Robert G.: Analytical Investigation of Fully Developed 

Laminar Flow in Tubes with Heat Transfer with Fluid Properties 
Variable Along the Radius. NACA TN 2410, 1951. 





10,000 


NACA TN 3145 


29 



0w 


= ,0 JO 


V^ c 

0,Q^/n _ 


0, 


0 , 


0, 


Figure 2. - Generalized temperature or concentration distributions for various Prandtl or Schmidt numbers. 


Nu = hD/k or Nu' = h'D/X 
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Figure 3. - Fully developed Nusselt numbers for heat or mass transfer against Reynolds 
number for various Prandtl or Schmidt numbers. 
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(a) Reynolds number, 10,000. 
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Figure 4. - Continued. Comparison of analytical and experimental results for fully developed heat and mass 
transfer. 
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Figure 4. - Concluded. Comparison of analytical and experimental results for fully developed heat and mass 
transfer. 







Figure 5. - Comparison of various analyses. Reynolds number, 10,000 
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velocity distribution. 
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Figure 6. - Continued. Variation of* local Nusselt number over fully developed Nusselt number for heat or mass transfer 
with X/D and Reynolds number. Uniform wall heat flux, uniform initial temperature distribution, and fully developed 
velocity distribution. 
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Figure 6. - Continued. Variation of local Nusselt number over fully developed Nusselt number for heat or mass transfer 
with X/D and Reynolds number. Uniform wall heat flux, uniform initial temperature distribution, and fully developed 
velocity distribution. 
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Figure 6. - Continued. Variation of local Nusselt number over fully developed Nusselt number for heat or mass transfer 
with X/D and Reynolds number. Uniform wall heat flux, uniform initial temperature distribution, and fully developed 

vpI nn-l fv * if 
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Figure 6. - Concluded. Variation of local Nusselt number over fully developed Nusselt number for heat or mass transfer 
with X/D and Reynolds number. Uniform wall heat flux, uniform initial temperature distribution, and fully developed 
velocity distribution. 
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Figure 7. - Continued. Variation of average Nusselt number over fully developed Nusselt number for heat or mass transfer 
with X/D and Reynolds number. Uniform wall heat flux, uniform initial temperature distribution, and fully developed 
velocity distribution. 
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Figure 8. - Concluded. Generalized distribution for liquid with variable viscosity. Fully developed flow; Prandtl 
ber, 10; \x/[i Q = (t/t 0 )~ 4 . (Vertical line is dividing line between eqs.) 









Figure 10. - Friction factors against Reynolds numbers for various values of heat-transfer parameter for liquids with 
variable viscosity. Fully developed flow; Prandtl number, 10; h/m-q = (t/tQ)’ 4 . 
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Figure 12. - Comparison of generalized distributions for gases with variable properties. Prandtl number, 0.73. (Vertical 
line is dividing line between eqs . ) 
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Figure 12. - Concluded. Comparison of generalized distributions for gases vith variable properties. Prandtl number 
0.73. (Vertical line is dividing line between eqs.) 


